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Abstract 

We introduce the Schlesinger transformations of the Gambier equation. The latter 
can be written, in both the continuous and discrete system of two coupled 

Riccati equations in cascade involving an integer parameter n. In the continuous case 
the parameter appears explicitly in the equation while in the discrete case it corresponds 
to the number of steps for singularity confinement. Two Schlesinger transformations are 
obtained relating the solutions for some value n to that corresponding to either n + 1 
or n + 2. 



1. Introduction 



The existence of Schlesinger transformations is one of the very special properties 
of Painleve equations [1]. These transformations are a particular kind of auto-Backlund 
transformations [2]. The latter relate a solution of a given equation to a solution of 
the same equation but corresponding to a different set of parameters. Schlesinger 
transformations do just that but the changes of parameters correspond to integer or 
half-integer shifts in the monodromy exponents. Both continuous and discrete (whether 
difference or multiplicative) Painleve equations have been shown to possess Schlesinger 
transformations [3]. For the discrete case (and in particular for g-Painleve's) the re- 
lation of Schlesinger transformations to monodromy exponents is not quite clear and 
their derivation requires both experience and intuition, in particular in the choice of 
the proper parameters. With this minor caveat we hasten to say that we do possess a 
systematic approach to the construction of Schlesinger transformations. It is based on 
the bilinear formalism [4] which can be used to construct Miura transformations [5] , the 
iteration of which can lead to auto-Backlund's. 

The logical conclusion of the above introduction is that Schlesinger transformations 
should exist only for Painleve equations. This is almost true with one exception. As 
we shall report in this paper there exists one equation which, without being a Painleve, 
does possess Schlesinger's. This equation is known under the name of Gambier (who 
first derived it) : it is the most general second-order ODE of linearisable type [6] . What 
makes possible the existence of Schlesinger transformations for this equation is the fact 
that its general expression involves an arbitrary integer n. It turns out that we can 
relate the solution of the equation for some value of n to that corresponding to n + 2. 
These transformations are, for the linearisable case, the analogues of the Schlesinger's. 
Moreover, as we shall show, the same procedure can be followed in the discrete case i.e. 
for the Gambier mapping [7,8]. 

In section 2 we shall review some basic facts about the Gambier equation. The 
Schlesinger transformations will be given in section 3 while section 4 is devoted to the 
study of the discrete case. 

2. The Gambier equation 

The Gambier equation is given as a system of two Riccati equations in cascade. 
This means that we start with a first Riccati for some variable y 

y' = -y 2 + by + C (2.1) 

and then couple its solution to a second Riccati by making the coefficients of the latter 
depend explicitly on y: 

x' = ax 2 + nxy + a. (2.2) 

The precise form of the coupling introduced in (2.2) is due to integrability require- 
ments. In fact, the application of singularity analysis shows that the Gambier system 



cannot be integrable unless the coefficient of the xy term in (2.2) is an integer n. This 
is not the only integr ability requirement. Depending on the value of n one can find 
constraints on the a, 6, c, a (where the latter is traditionnally taken to be constant 1 
or 0) which are necessary for integrability. 

The common lore [9] is that out of the functions a, 6, c two are free. This turns 
out not to be the case. The reason for this is that the system (2.1-2) is not exactly 
canonical i.e. we have not used all possible transformations in order to reduce its form. 
We introduce a change of independent variable from t to T through dt — gdT where g is 
given by --^ = b^z^, a gauge through x = gX and also Y = gy — ^^f • The net result 
is that system (2.1-2) reduces to one where 6 = while a remains equal to or 1. It 
is clear from the equations above that n must be different from 2. On the other hand 
when n = 2 the integrability condition, if a = 1, is precisely 6 = 0. So we can always 
take 6 = 0. (As a matter of fact in the case a = an additional gauge freedom allows 
us to take both 6 and c to zero for all n, even for n = 2). Thus the Gambier system can 
be written in full generality 

y' = -y 2 + c (2.3a) 

x = ax 2 + nxy + a. (2.36) 

One further remark is in order here. The system (2.3) retains its form under the 
transformation x — > 1/x. In this case n — > — n and a and —a are exchanged. Thus in 
some cases it will be interesting to consider a Gambier system where a is not constant 
but rather a function of t. Still, it is possible to show that we can always reduce this 
case to one where a = 1, while preserving the form of (2.3a) i.e. 6 = 0. To this end 
we introduce the change of variables dt = hdT, x = gX and Y = hy — with 
h = a 2 /( n ~ 2 \ g = a n /( n ~ 2 \ With these transformations system (2.3) reduces to one 
with a = 1 and 6 = 0. (In the special case n = 2, with 6 = integrability implies 
a =constant, whereupon its value can always be reduced to 1). 

3. SCHLESINGER TRANSFORMATIONS FOR THE GAMBIER EQUATIONS 

The theory of auto-Backlund transformations of Painleve equations is well estab- 
lished. As was shown in [2] the general form of auto-Backlund transformations for most 
Painleve equations is of the form: 

ax' + fix 2 + 7X + 5 . . 

x = — ; — tt^ w I 3 - 1 ) 

ex' + (x 2 + i]x + 9 

In the case of the Gambier equation considered as a coupled system of two Riccati's it 
is more convenient to look for an auto-Backlund of the form: 

_ _ axy + f3x + >yy + 8 

X ~ (Cy + V )(0x + K) ■ [6 - Z) 



with a factorized denominator, with hindsight from the discrete case. We require that 
the equation satisfied by x do not comprise terms nonlinear in y. We examine first the 



case ( ^ and reach easily the conclusion that there exists no solution. So we take 
C = 0, i] = 1 which implies that a and 7 do not both vanish (otherwise (3.2) would have 
been independent of y). We find in this case a = and thus the general form of the 
auto-Backlund can be written as: 

From (3.3) we can obtain the two possible forms of the Gambier system auto-Backlund: 

x = [3x + <yy + 5 (3.4) 

x = ^±J1±1. (3.5) 
x + k J 

As we shall see in what follows both forms lead to Schlesinger transformations. 

Let us first work with form (3.4). Our approach is straightforward. We assume 

(3.4) and require that x, satisfy an equation of the form (2.3b) while y is always the 

same solution of (2.3a). The calculation is easily performed leading to: 

a = 72H — r^H , (3-6) 

Ti + l n 

where 7 satisfies: 

^ = (3-7) 

7 n + 2 a v ' 

Here we have assumed a ^ otherwise x, does not depend on x and (3.6) does not 
define a Schlesinger. The parameters of the equation satisfied by x are given (in obvious 
notations) by: 

ft + n + 2 = (3.8a) 

a = ^ (3.86) 

7 

and 

/ aa la" n + 3 a' 2 \ /0 n . 

(7 = 7 c+ - + t 3.8c) 

V n + 1 n + 2 a (n + 2) 2 a 2 / 

Thus (3.6) is indeed a Schlesinger transformation since it takes us from a Gambier 
system with parameter n to one with parameter n = — n — 2. It suffices now to invert x 
in order to obtain an equation with parameter iV = n + 2. Expressions (3.6) and (3.8) 
can be written in a more symmetric way: 

a' 



na 

and 

n + 1 

~d a' 



ax — ax = (n + l)(y — — ) (3-9) 

n+1 = -(n+1) 
„a' a' 

n ^a =n ~a . (3-10) 

/ t \ f 1 / a \ ' 1 a' 2 \ 
7 =( n+ l)( c __(_) +— — . 

V n v a ' n 2 a z / 



The inverse transformation can be easily obtained if we introduce 7 such that 
07 = — (n + 1) = —07. We thus find 

x = y*{+ -J^-x+Z- (3.11) 
n + 1 n 

and the relations (3.10) are still valid. 

Iterating the Schlesinger transformations one can construct the integrable Gambier 
systems for higher n's and obtain by construction the functions which appear in them. 
However it may happen that when we implement the Schlesinger we find a = 0. If we 
invert x we get a system with N = — h = n + 2 but A = for which one cannot iterate 
the Schlesinger further. 

Let us give example of the application of this Schlesinger transformation. Let us 
start from n = 0, in which case we find n = —2 and, after inversion, N = 2. For n = 
we start from a = — 1 and a = or 1 (always possible through the appropriate changes 
of variable). This leads to a = —1, a = —c + a and the Schlesinger reads: x = —y + x. 
Next we invert x and have X = 1/ (x — y) . We find thus that the Schlesinger takes us 
from 

y' = -y 2 + c 

x' = -x 2 + a 

to the system 



y' = — y 2 -\- c 

X' = AX 2 + 2Xy + E (3 - 13) 
with A = c — a, E = 1. In the particular case n = 2, a change of variables exists which 
allows us to put A = — 1 (unless A = 0), without introducing b in the equation for y, 
while keeping E = 1 and changing only the value of c. Thus the generic case of the 
Gambier equation for n = 2 can be written with A = — 1. Eliminating y between the 
two equations we find: 

/2 3 1 

x " = — - 2xx' - - + (2c+l)x. (3.14) 

2x 2 2x K 1 V 1 

This is the generic form of the n = 2 Gambier equation and it contains just one free func- 
tion. The nongeneric cases corresponding to A = and a = or 1 can be constructed 
in an analogous way. 

We now turn to the second Schlesinger transformation corresponding to the form 
(3.5). As we shall show, a Schlesinger transformation of this form does indeed exist and 
corresponds to changes in n with An = 1. Let us start from the basic equations (2.3). 
Next we ask that x defined by (3.5) indeed satisfy a system like (2.3). We find thus 
that and k = —x$ and 7 must be given by: 

— = Vo + ——r 3.15) 
7 n + 1 

where yo is a solution of the Riccati (2.3a) and xq a solution of (2.3b), obtained with y 

replaced by yo. We introduce the quantities xq = a = — 2 ^ a . In this case (3.15) 
becomes: 

i . 2ax 2x x . . 

— = y + = y H , 3.16 

7 n+1 7 



where 

n + n + 1 = 0. (3.17) 
We have thus, starting from a generic solution x, y of (2.3) for some n, the Schlesinger: 

x = x + 1{V ~ VQ) (3.18) 

X — Xq 

where x is indeed a solution of (2.3) for n = — n — 1 for the same y 

x = —ax 2 + nxy + a (3.19) 



where a has been defined as — 2£o anc [ 

7 



ry / _L_ \£ \ 

a = —!— (a' + a 2 x + ay (n + 2) . (3.20) 

n + 1 V n + 1 / 

Note that is a solution of the same equation with y replaced by yo- As in the previous 
case if we invert x we obtain an equation corresponding to iV = n + 1. 

As an application of the An = 1 Schlesinger we are going to construct the n = 1 
equation starting from the n = case, i.e. system (2.3) with n = 0. From (3.18), the 
Schlesinger reads 

5 = 7 (a+^l^) (3.21) 

X — Xq' 

where 7 satisfies the differential equation 

l' 

— = y + 2ax (3.22) 

7 

and where y is a particular of (2.3b) and x is a solution of (2.3a) with n = 0. Then x 
satisfies a Gambier equation in which n = —1, a = 7(0' + 2a 2 xo + 2ayo) and a = 0. To 
obtain the n = 1 Gambier equation, we define X = 1/x and we arrive at the following 
system 

y' = -y 2 + c 

(3.23) 

X' = AX 2 + Xy 

where A = —a. The system (3.23) is the generic n = 1 case since, in this case, the 
condition for (2.3) to be integrable is a = 0. 

It is worth pointing out here that the Schlesinger transformation corresponding to 
An = 2 was known to Gambier himself. As a matter of fact when faced with the problem 
of determining the functions appearing in his equation so as to satisfy the integrability 
requirement, Gambier proposed a recursive method which is essentially the Schlesinger 
An = 2. On the other hand the Schlesinger An = 1 is quite new and we have first 
discovered it in the discrete case whereupon we looked for (and found) its continuous 
analogue. 



4. The Gambier mapping 



The Gambier equation has been examined already in [7,8] and its discrete equiva- 
lent has been proposed there. These constructions of the Gambier mapping were ad hoc 
ones, in the sense that we assumed a form and implemented the singularity confinement 
discrete integrability criterion in order to obtain the integrability conditions. In what 
follows, we shall use a slightly different approach based on the singularity structure. 

Our starting point is the discrete equivalent of the system (2.3). We have thus 
one equation which is the discrete analogue of the Riccati, i.e. a homographic mapping 
for y and another homographic mapping for x, the coefficients of which depend linearly 
on y. Our derivation will be based on the study of singularities of the system. The 
general homographic equation for y involves three free parameters, but since we have 
the freedom of choice of a homographic transformation on y, we can always reduce it 
to y =constant (i.e. y = y). However, in the system under study, our aim is to study 
the singularities of x induced by special values of y. One could choose the singularity to 
enter at point no if the value of y has some special value depending on no, say, /(no). 
This would introduce one function in the homographic mapping. However what is even 
more convenient is to decide what the special value of y is, say y = for all n, at the 
price of the loss of part of the homographic freedom. Then the special value will 
occur for some n depending on the particular solution. Of course, if we allow the full 
homographic freedom we are back to the starting point i.e. with three free functions. 
However we decide to have only one free function and thus we simplify the mapping by 
chosing its form so that it presents the pattern {—1, 0, 1}. This fixes two of the functions 
and the result is: 

where c is a function of n and we use the notations y = y(n), y = y(n + 1). 

Next, we turn to the equation for x. This equation is homographic in x. However 
we require that when y takes the value 0, the resulting value of x be oo. Thus the 
denominator must be proportional to y, and since we can freely translate x we can 
reduce its form to just xy. The remaining overall gauge factor is chosen so as to put 
the coefficient of xy of the numerator to unity resulting to the following mapping: 

x = x(y-r) + q(y-s)^ 
xy 

The system (4.1-2) is a discrete form of the Gambier system. In order to study the 
confinement of the singularity induced by y = we introduce the auxiliary quantity i/jn 
which is the iV'th iterate of y = in equation (4.1), N times downshifted. Thus ipo = 0, 

c+c 

ipi = c, ip2 = ^pj, etc... The confinement requirement is that after N steps x becomes 
in such a way as to lead to 0/0 at the next step. Thus the mapping (4.2) has in fact 
the form: 

_ = x(y-r)+ q (y-^ N )_ 
xy 



Thus when at some step N we have y = V>jv and x = 0. On the view of (4.3) x will 
then be indeterminate of the form 0/0. However it turns out that, in fact, this value is 
well-determined and finite. Let us take a closer look at the conditions for confinement. 
The generic patterns for x and y are: 

y:{ Vi _ V>2 • • • V'iv } 

x : { free oo —. — ... free }. 

At N = 1 it is clearly impossible to confine with a form (4.3) since we do not have 
enough steps. In this case the only integrable form of the x-equation is a linear one. 
The first genuinely confining case of the form (4.3) is N = 2. From the requirement 
x = we have r = ifti and q free: this is indeed the only integrability condition. For 
higher iV's we can similarly obtain the confinement condition which takes the form of 
an equation for r in terms of q. 

At this point it is natural to ask whether the mapping (4.1)-(4.3) does indeed 
correspond to the Gambier equation (2.3). In order to do this we construct its continuous 
limit. We first introduce: 

c = e 2 D 

eD ( 4 - 4 ) 

y 



Y + H 

with H a* D'/(2D) and obtain the continuous limit of (4.1) for e — > 0. We find as 
expected 

Y' = -Y 2 + C (4.5) 

i.e. eq. (2.3a), where C = D — ^ + \%r- Using (4.4) and (4.1) we can also compute 
ipN and we find at lowest order: 

iJj N =e 2 y N with ^ N w N(D-e^-D') +e 2 $ N . (4.6) 

where is an explicit function of D depending on N . 
Next we turn to the equation for x and introduce: 

r = e 2 R 

1 e RD' 
X ^2 + 2X-^ (4.7) 

1 9 

q w -- + e 2 Q 

and for the continuous limit of the form (2.3b) to exist in canonical form (i.e. 6=0, 
<7=1) we find that we must have 



(4.8) 



This leads to the equation for x: 



X' = AX 2 + NXY + 1 (4.9) 

with A = + 1)^2 — 4Q. Moreover the confinement constraint implies 

a differential relation between D and Q which depends on N. We can verify explicitly 
in the first few cases that this is indeed the integrability constraint obtained in the 
continuous case. For instance, for N = 2, just imposing (4.8) in order to have the 
canonical form 6=0, cr=l, is sufficient for integrability. 

Once the singularity pattern of the Gambier mapping is established we can use it in 
order to construct the Schlesinger transformation. Let us first look for a transformation 
that corresponds to AN = 2. The idea is that given the iV-steps singularity pattern 
of the equation for x we introduce a variable w with N + 2 singularity steps where we 
enter the singularity one step before x and exit it one step later. The general form of 
the Schlesinger transformation, which defines w, is: 

w = X y -^ N+ \ (4.10) 

y 

where X is homographic in x. The presence of the y and y — i^>n+i terms is clear: they 
ensure that w becomes infinite one step before x, and vanishes one step after x. Next 
we turn to the determination of X. Since X is homographic in x we can rewrite (4.10) 
as: 

w _ ^x + (3y-jj N+1 ^ 
y "fx + 5 

Our requirement is that w becomes infinite when y = for every value of x. This 
statement must be qualified. The numerator ax + (3 will vanish for some x (namely 
x = —(3 /a) so this value of x must be the only one which should not occur in the 
confined singularity. Indeed there is a unique value of x where instead of being confined, 
the singularity extends to infinity in both directions of the independent variable n, while 
the only nonsingular values of the dependent variable occur in a finite range. The value 
of x such that x is finite and free even though y is zero is such that the numerator 
—xr — qtf)N of x vanishes. For this value of x, the values of the dependent variable are 
fixed for n < and n > N + 1 and the value can be considered as 'forbidden'. Thus 
ax + [3 = xr + qipN up to a multiplicative constant. Similarly when y = V'jv+i, w must 
vanish. Thus 7^ + 5 must not be zero except for the unique value of x that does not 
occur in the confined singularity. Note that y = i^n+i means y = ip and the only 
value of x that comes from a nonzero x in that case is x = (ip N — r)/i/j . In that case 
the values of the dependent variable are fixed for n > and n < —N — 1. This value of 
x being 'forbidden', -fx + 5 must be proportional to ip N x — — r). We now have the 
first form of the Schlesinger: 

w _ xr + qjj N yj-^N+i 

y ± N x + l-± n 



where the proportionality constant has been taken equal to 1 (but any other value 
would have been equally acceptable). Here w effectively depends on x unless r(r — 
tp ) = qip ipN- But in this case the mapping (4.3) is in fact linear in the variable 
£ = (x — 1 + r/V'jv) -1 - This case is the analog of the case a = in the continuous case 
where the Schlesinger does not exist. 

Let us give an application of the Schlesinger transformation by obtaining the N = 2 
equation starting from N = 0. We have always the equation for y which reads: 

and ip = 0, V'l = c. For N = the equation for x reads: 

__ x(y-r)+qy _ x + q 
xy x 

since for integrability r = and indeed N = means that the x equation does not 
depend on y. We introduce the Schlesinger: 

w = x V -^ (4.15) 

y 

(where we first write (4.12) for arbitrary r and since r factors we take the limit r — > 
afterwards). Using (4.14) and (4.15) to eliminate x we obtain the equation for w: 

- n ,yw + q(y-c) 

w=(l — c) — -, r . (4.16) 

This equation is of the form (4.3) but not quite canonical. We can transform it to 
canonical form simply by introducing y instead of y because indeed w is infinite one 
step before x, so w = oo means x = oo i.e. y = 0. We obtain thus: 

__ w(y- c) +q(l + c)(y-^ 2 ) 
yw 

with ip 2 = (c + c)/(l + c) which coupled to (4.13) is indeed & N = 2 Gambier mapping. 

As we pointed out above the N = 1 case is not included in the parametrisation 
(4.1-4.3): the x-mapping must be linear in order to ensure integrability. Thus we are led 
to study the linear case separately. For an arbitrary N, the general form of the linear 
x-mapping can be obtained using confinement arguments in a way similar to what we 
did for the generic, nonlinear, case. We obtain 

_ = xiy^M+g 

y 

where g is free. The Schlesinger transformation is again given by: 

w = X V ' ijN+1 (4.19) 

y 



and arguments similar to those of the nonlinear case allow us to determine the form of 
the homographic object X leading to: 



M= ^-gy-^ + i (420) 
y x± N - g 

Thus one can perform a Schlesinger in the linear case. This is not in disagrement with 
the continuous case. It is, in fact, the analog of the case where a = but a ^ (which 
is linear in 1/x) for which the Schlesinger can be performed. The analog of the case 
a = and a = is the situation when g = kipN with constant k in which case the 
mapping rewrites £ = £(y — iPn)/u with £ = x — k. Then w does not depend on £ (or 
x) and (4.20) does not define a Schlesinger in analogy to the case r(r — ip ) = qifj N i^N 
in the nonlinear case. 

Using this form of the Schlesinger transformation we can, for example, construct 
the N = 3 case starting from the N = 1 case. In the case N = 1, the mapping for x is 
given by: 

x=(y-^ x + 9. (4.21) 

y 

Using equation (4.20), we introduce the Schlesinger: 

In order to simplify the final expression, we define p with g = ep. We then have the 
following equation for w: 

c(c -l)(cw ((p - p)y + c(p - p)) + c(p - p) (y(c + 1) - (c + c)) ) 

W= ^ r r J ~. (4.23) 

(p - p)cgw{y + c) 

To give this equation in the same parametrisation as (4.3), we first write it in terms of 
y and then use a multiplicative gauge u> = 4>w to put to unity the coefficient of uy in 
the numerator of u. We then have: 

m= "(v- r ) + (4.24) 
uy 

where _ 

(p - p) (p - p) (2g + c + 1) 

<2 = 



(c(p-p) +p-p)(g(p-£) +P-^) 

r = + (4.25) 

c(p — p) + p — p 
— CC + c + c + c 

^ 3 ~ 2g + c + 1 

Finally, we examine the possibility of the existence of a AiV = 1 Schlesinger. In 
this case, the structure of the transformation will be obtained by asking that the N + 1 



case enter the singularity one step before the N case but exit at the same point. The 
general structure is thus: 

rx + qi(j N y-r] 
w = (4.26) 

where rj and £ must be determined. We do this by requiring that the equation for 
w contain no coefficients nonlinear in y. As a result we find that r\ must satisfy the 
equation (4.1) for y: 

V + c . 

5?= < 4 ' 27 > 

and £ the equation (4.3) for x with r\ instead of y: 

£ = g^Z ^ + ^ ~ M . (4.28) 

We remark here the perfect parallel to the continuous case (and as we pointed out the 
discrete case led the investigation back to the continuous one). Let us point out here 
that the w obtained through (4.26) does not lead to w = at the exit of the singularity 
(i.e. when x = 0, y = tpx) and a translation is needed. One has in principle to define a 
new variable 

q 

u = w - w(x = 0, y = tpN) = w + -(V'jv - rj). 

We are now going to study the particular case where we construct N = 1 starting 
from N = 0. As in the continuous case our starting point is the = nonlinear 
equation. Thus we are starting with the decoupled equation (4.14) for x. We write the 
Schlesinger as: 

W = V ~ T1 X t . (4.29) 

Once more we find that r\ must satisfy (4.27) while £ satisfies £ = 1 + q/£ i.e. the same 
equation as x. Using (4.27) we can easily obtain the equation for w corresponding to 
N = 1. We write this equation in terms of y in order to enter the singularity with y = 0. 
The expression of w is indeed linear in w and reads 

_ (g + Qw(c -y)+ q(y(r] + 1) - c - rj) 

w = — ; — . (4.3UJ 

yq( V + l) 

In order to cast it in the parametrisation of (4.21), we introduce a shift in the w: 

ui = w + ujo. (4-31) 

We require that in the numerator of the equation for u, y appears only as a product 
with u>. We find that ujq must satisfy 

= _Mi + + «(■>+!) (4 32) 

q{V+l) 



and we obtain the following equation for u 



_ = (g + CMg - y) - c^ (q + Q - q(c + 77) 33 
2/9(77 + 1) 

which is in the form of (4.21) up to a multiplicative factor in to. We note that ujq plays 
the role of x in equation (3.18) in the continuous case. Indeed, ujq does satisfy (4.33) 
for y = 77. 

Finally we derive the AN = 1 Schlesinger for the case of a linear mapping (4.18). 
We start from: 

w = (4.34) 

y x-t, 

and again require for w an equation with coefficients linear in y. We find that 77 must 
again be a solution of the equation for y i.e. it must satisfy (4.27) and moreover £ is a 
solution of (4.18) with y = 77: 

{= fcM + g . (4.35) 
Thus the list of the Schlesinger transformations of the Gambier mapping is complete. 
5. Conclusion 

In this paper we have shown that the Gambier system possesses Schlesinger trans- 
formations just like the Painleve equations. This is a most interesting result given that 
the Gambier system is C-integrable (in the terminology of Calogero [10]) i.e. integrable 
through linearisation, and not S-integrable. 

We discovered that both the continuous and the discrete systems possess two 
kinds of Schlesinger transformations: one that allows changes of N by two units and 
one where the changes of N are by one unit. In the discrete case our approach was 
based entirely on the singularity confinement aproach. We have shown that a study of 
the singularities allows one to determine the form of the Gambier mapping and at the 
same time its Schlesinger transformations. This is one more argument in favour of the 
singularity analysis approach to the study of discrete systems. 
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